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Abstract 

We consider the model of "Chain Inflation," in which the period of inflation in our 
universe took the form of a long sequence of quantum tunneling events. We find that in 
the simplest such scenario, in which the tunneling processes are uniform, approximately 
10^ vacua per e-folding of inflation are required in order that the density perturbations 
produced are of an acceptable size. We arrive at this conclusion through a combina- 
tion of analytic and numerical techniques, which could also serve as starting points for 
calculations with more general sets of assumptions. 



1 Introduction 



In 1980, Guth proposed the inflationary universe scenario to solve the horizon and flatness 
problems [1]. In this picture, the very early universe experienced a period of exponential 
expansion driven by the vacuum energy of a scalar field. In the original model, the scalar field 
was supposed to be stuck in a false vacuum of its potential. After some time, however, quantum 
mechanical tunneling would bring the field to the true minimum of the potential and inflation 
would come to an end. 

The tunneling process in quantum field theory was originally described by Coleman [2]. A 
tunneling event is characterized by the formation of a "bubble," in which the inside of the 
bubble is located in the new true-vacuum phase, while the outside of the bubble is located 
in the old false-vacuum phase. The vacuum energy lost in the interior of the bubble due to 
its formation is stored in the bubble wall; energy is conserved in this process. The bubble 
subsequently expands, at a speed quickly approaching the speed of light, eating up the false 
vacuum in favor of the true one. When two bubbles collide the energy stored in the walls is 
converted into a spray of particles, which eventually thermalizes. 

Unfortunately, Guth's original model had a fatal flaw: If the tunneling rate from the false 
vacuum was too large, then there was too short a period of inflation to yield the 60 e-foldings 
of expansion necessary to solve the horizon and flatness problems. If the tunneling rate was 
too small, then the true vacuum bubbles formed too far away from each other to ever have 
consistent collisions and percolation. The interiors of the bubbles would then tend to remain 
forever empty, with the energy being stored permanently in the bubble wallsjl] It was argued in 
[3] that there is likely no intermediate value for the tunneling rate which can yield an acceptable 
number of e-foldings while also allowing bubbles to collide and percolate. 

As a result of this problem, models of inflation based on slowly rolling fields [1] [5], rather 
than vacuum tunneling, have become the standard paradigm. On the other hand, Freese and 
Spolyar, in their paper "Chain Inflation" [B], proposed a simple way to extend and save Guth's 
false vacuum style inflation. Their idea was to choose a tunneling rate on the high side, so 
that bubbles would indeed percolate, but at the same time have a long sequence, or "chain," of 
vacua to ensure an acceptable number of e-foldings. In each step of the chain, only a fraction 
of an e-folding is obtained, but overall, inflation can be sustained for a sufficiently long time. 
The biggest open question in this model is whether or not it can lead to a phenomenologically 

^Bubble collisions do occasionally happen in this scenario, but the inhomogencities in the resulting radiation 
tend to be much too large. 
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acceptable spectrum of density perturbations. The density perturbations need to be scale 
invariant, of an appropriate amplitude, and have very small departures from gaussianity. It is 
not at all clear how such a spectrum could arise in chain inflation. 

As noted in [6] the main source of density perturbations in chain inflation should come from 
the intrinsic randomness of the tunneling process itself. Different points in space will reach the 
final step in the chain at different times, leading to inhomogeneities. Scale invariance can most 
easily be satisfied by assuming a time translation symmetry amongst the tunneling events - i.e., 
by assuming all steps in the chain are essentially identical. An example model with this feature 
is an axion with a tilted cosine potential, as discussed in [7] . In this paper we will concentrate on 
this case, although we will have comments concerning extensions to non-uniform tunneling as 
well. Calculating the density perturbations analytically in chain inflation is extremely difficult, 
and we will instead use computer simulations for much of our analysis. 

During the period of inflation responsible for the observed modes of the density pertur- 
bations, the Hubble parameter H may be taken as roughly constant. Then the only free 
dimensionless parameter on which the amplitude of the perturbations can depend is r~^if^, 
where F is the tunneling rate per unit volume for the uniform tunneling events. The constraint 
of bubble percolation requires that T-^H^<1, but since the amplitude of the perturbations 
is required to be ~ 10~^ by the COBE normalization [8], naively we might expect to find a 
requirement of T~^H'^ <^ 1 instead. This is equivalent to having many vacuum transitions per 
e-folding of inflation. Indeed, we will ultimately show that about 10^ transitions per e-folding 
are needed in order to obtain sufficiently small density perturbations. 

In section 2 we describe our specific assumptions in more detail, as well as outline the 
method by which the density perturbations will be calculated. This includes a demonstration 
that the amplitude of perturbations will be given in terms of the RMS fluctuations of the 
vacuum number at a point. In section 3 we give a description of the algorithm used in our 
simulations. In section 4 we give our results, and we summarize and consider future directions 
in section 5. Further details concerning the simulations and data analysis can be found in the 
appendices. 

2 Density Perturbations 

The generation of density perturbations in chain inflation is not entirely unlike the situation 
in ordinary slow-roll inflation. In both cases, the perturbation modes can be understood as 
arising from x-dependent field fluctuations, generated by the quantum mechanical randomness 
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inherent in the microscopic evolution of each system. Unless stated otherwise, we shall work 
in synchronous gauge, in which "t" labels proper time, and "x" labels worldlines of free-falling 
observers, with the physical distance between two points being given at zeroth order in the 
perturbations by a|x — y|, where a is the scale factor, "t = 0" may be taken to denote a 
spatially flat hypersurface in the distant past at a time at which our present horizon scale 
was located within a single bubble. Note that we shall ignore the back-reaction of the metric 
perturbations on the evolution of the infiaton field, as this would be a higher order effect. 

In the slowly rolling case, a given mode starts out far in the past in a coherently-evolving 
quantum ensemble - a component of the scalar field vacuum in de Sitter space. As the mode 
inflates outside of the horizon, the ensemble decoheres, leaving distinct "classical" field fluctu- 
ations |9] [To] [H]. The field then rolls down its potential deterministically throughout space, 
but its fluctuations cause nonzero relative timings 



where (j)o{t) is the spatially averaged field. These in turn lead to metric and matter density fluc- 
tuations. As observers looking at the cosmos today, we can then understand the fluctuations in 
the microwave sky and matter distribution as ultimately arising from the quantum-mechanical 
uncertainty of field values in the de Sitter vacuum. 

Similar considerations also apply in chain inflation, although the dynamics are significantly 
different. As the tunneling field proceeds down the chain and the universe undergoes successive 
periods of percolation, the inherently random placement of bubble nucleation sites across space- 
time causes certain patches of space to advance further along the chain than others, translating 
into timing fluctuations. As in the slow-roll case, these timing fluctuations are produced on 
scales smaller than the Hubble radius, and evolution on trans- Hubble scales is again essentially 
deterministic. The density perturbations thus produced can again be seen to be a product of 
quantum randomness, but now originating in the tunneling process instead of the decoherence 
of de Sitter vacuum modes. 

Unfortunately, finding the statistical properties of the evolution of a given bubble network 
can be highly nontrivial. The evolution is tied to the details of bubble wall collisions and the 
interaction with radiation from previous percolations. If the tunneling rate varies wildly from 
vacuum to vacuum (as might be expected due to the exponential sensitivity of the tunneling 
rate to the form of the potential) then things can get very messy. On the other hand it is 
possible to make a well defined set of simplifying assumptions so that a calculation of density 
perturbations becomes tractable. We will assume the following four conditions: 




0o(t) 
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1. The size of the bubble immediately after tunneling is small compared to its characteristic 
tunneling time, r = T~^^^. 

2. Upon collision, bubble walls dissipate quickly into radiation, on time scales shorter than 
(or at least not much longer than) r. 

3. Tunneling which skips steps in the chain is sufficiently suppressed that it can be ignored. 

4. All steps in the chain have identical tunneling events. This includes tunneling rates (per- 
haps effective thermal tunneling rates), as well as the total vacuum energy lost at each 
step. 

As noted in the introduction, this final assumption is the only obvious way to ensure that 
the final spectrum is scale invariant, and it also greatly simplifies our calculation method. To 
the extent that our assumptions are not realistic for a given scenario, our results may be viewed 
as a baseline for more complicated treatments and/or construction of more advanced simulation 
methods. 

Note that as the universe expands, radiation left over from earlier percolations is infiated 
away. However, it takes about an e-folding of infiation for a given quantity of radiation to be 
noticeably diluted. This implies that there is always present a bath of radiation, which has 
an average energy density given by the total vacuum energy lost in about one Hubble time. 
Eventually the vacuum energy drops below the energy density in this radiation, and infiation 
comes to an end. Note that assumptions 1 through 4 need only be valid up to an epoch about 
30 to 60 e-foldings before the end of infiation in order to satisfy the observational constraints 
from the CMB (i.e., up to the epoch when our current Hubble scale exited the horizon). If they 
are always true, on the other hand, then in fact we can make an estimate for the reheating 
temperature at the end of infiation: In this case the Hubble parameter may be written as 



where pr is the average radiation energy density, and pv is the vacuum energy. Since we have 
about one e-folding worth of radiation at any given time, pr satisfies 

~ (3) 

with e representing the vacuum energy lost in each step of the chain|§ This implies the reheating 
temperature will be about 



^The naive expectation that the vacuum energy goes like pv PV:i ~ ^t/r will be confirmed later in our 
simulations. 
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Finally, we note that in what follows, we will treat the Hubble parameter as roughly constant 
during the period of inflation responsible for generating the modes of the density perturbations 
probed by observations of the CMB. This should be a good approximation, since this period 
is expected to be fairly short compared to the epoch of inflation remaining after our present 
Hubble volume exited the horizon. This implies that during the era of production of the 
observed density perturbations, the vacuum energy could only have changed by a small fraction. 

With our setup and assumptions now in hand, we may next lay out a systematic program 
for calculating the size of the density perturbations in this model. In order to preserve the 
analogy with slow roll inflation, we will denote the vacuum number as a function of position 
by (/)(x, t). Note, however, that with this definition is a dimensionless field. 

We begin by constructing the spatial two point correlator of fluctuations, defined as the 
expectation value (((^(x, t) — {4>{t))) {(j){y,t) — {4>{t)))) = (0(x, t) 0(y, t)) — (</)(t))^, where we 
average over an ensemble of universes (or physically remote regions in a single universe). Note 
that the zero-mode of this quantity is unobservable, as it will incorporate fluctuations of (f) 
inaccessible to any single observer viewing only a particular local patch. Other Fourier modes 
of this object do have a well-defined physical meaning, however; they represent the power in 
each field fluctuation mode, and will ultimately translate into the power of density fluctuations 
observed in the CMB. 

We will now demonstrate that the two point correlator of (p fluctuations may be expressed 
on large scales in terms of the standard deviation of at a single point, plus a constant term 
which contributes only to the zero-mode. In order to see this, let us consider the evolution of 
the correlator along two worldlines at positions x^^ and x^ which are separated by a physical 
trans-Hubble distance L at the time t. There will be an earlier time 



before which the two worldlines were closer together than the characteristic size of the percolat- 
ing bubbles, and therefore typically inside the same bubble. Let us call the common field value 
at this time (j){tr{L)). As the worldlines are driven apart by the Hubble flow, they depart from 
their mutual bubble and experience distinct but correlated evolutions until a time tu shortly 
after they leave each other's horizons. Their evolution subsequently decouples. The complete 
evolution therefore breaks up into three distinct periods: 



tr{L) r^t-H-^ log 



(5) 
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1. t < 1^.: Completely correlated evolution. At the end of this period, worldlines A and B 
share some common value (j){tr{L)). 

2. t-r < t < tn'- Partially correlated evolution. Each worldline acquires a correlated advance 
in field value and 6b, respectively, with (SaSb) — (5)^ > 0. 

3. t > tH'- Completely uncorrelated evolution. Each worldline acquires an uncorrelated ad- 
vance in field value and A^, respectively, with (A^ A^) — (A)^ ~ and {Aa/b ^a/b) — 
(A) (5). 

The situation is illustrated in figure [TJ Representing ^^(t) = (f){tr{L)) + 6a + A^, and 
similarly for B, we obtain 

{Mt) Mt)) - m)? = [mAL)?) - m^{LW] + [{Sa sb) - [sf] . (6) 

If we consider a different value of L at the same t, the argument repeats. However, we are 
assuming that the system has nearly-constant r and nearly-constant H, and that its evolution 
is essentially insensitive to absolute time. In particular, period #2 above - evolution from r 
separation to Hubble separation - is independent of L, and so should look essentially identical 
for all pairs of points that end up outside the horizon by the time t. The generic trans-horizon 
spatial correlator then takes the form 

(Mt) Mt)) - m)? = MtriL))') - m^{LW] + const, (7) 

with the constant referring to terms independent of L. As claimed, the problem of deter- 
mining the two-point correlator of thus translates into the problem of determining its RMS 
fluctuations at a point function of time. 

If the system is truly scale-invariant, then different slices of time evolution should behave 
statistically in the same „av. accruing identical fluctuationsi This implies that {4>{t)) and the 
RMS fluctuations of (p take the form 

= a{Hr) Ht, (8) 

ct^RMsitf = mf) - m)? = PiHr) Ht, (9) 

^There will be very early times before which the scale of our entire observable universe was contained within 
a distance scale less than t, and during which the Hubble parameter could have been appreciably higher than 
its later size. Such periods will contribute to (pRMsit)'^ in ways not proportional to t, but their effect will be 
L-indcpendent on scales within our present horizon. 
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Figure 1: Schematic representation of the evolution of the 2-point correlator of fluctuations. 
(pitr) represents the totally correlated evolution at scales less than r, 5's represent the partially 
correlated evolution on scales between r and 1/H and A's represent the totally uncorrelated 
evolution on scales greater than 1/H. 

with a and /3 yet-to-be-determined functions of the nucleation density. Extra factors of H have 
been inserted for convenience of normalization. It now follows from equations [5], [7] and [9] that 



We repeatedly verify the behavior [H] and M hi our numerical simulations, and explicitly demon- 
strate the logarithmic scaling of the correlator for several values of r and t. We describe the 
details of these simulations in the subsequent sections. Now, taking the Fourier transform of 
the correlator (and invoking spatial homogeneity) we obtain 



(0(x, t) 0(y, t)) - {<f){t)f = -(3{Ht) log |x - y I + const. 



(10) 



(50k(t)50k'(^)) = (27r)=^(53(k-k') 




(11) 



which is the desired scale-invariant form (here the k's label comoving wavenumbers). This 
yields a power spectrum for 5(j) given by @ 

V5^ = (3{Ht). (12) 

Now, when a mode of the density perturbations ^ re-enters the horizon during matter 
domination, its size on a comoving hypersurface (one on which the momentum density vanishes) 
is related to the spatial curvature on that surface, through the relation [12] 



^ = In,. (13) 

P 5 



1 ^(3) 



Here T^k is defined by 

(14) 

The observations of the CMB by COBE then constrain the current power spectrum V-ji, for 
modes of order the Hubble radius today. The value obtained is [8] 

{V^o^'^y^ =A.^xlQ~\ (15) 



This quantity may be related to the synchronous gauge Vsct, which appears in equation 
evaluated at horizon exit. To see this, consider the epoch when a perturbation mode first exits 
the horizon during inflation. In synchronous gauge, the fluctuations in the vacuum number yield 
fluctuations in the spatial curvature through Einstein's equations of order TZq ~ ~ ~^p$' 
with the vacuum energy py taking the form pv = pv,i—^4>- We then must make a transformation 
to comoving coordinates and see how the spatial curvature changes. 

In slow roll inflation, comoving hypersurfaces correspond on trans-horizon scales to surfaces 
of constant field value [12], and the same will be true in chain inflation. The momentum density 
of a chain-inflating field is localized inside bubble walls and is constantly downconverting into 
radiation. It is therefore the radiation which will dominate the momentum density of the 
system, and we must choose a frame in which that momentum averages out to zero. This 
will happen in a spatial slicing where tunneling to a particular vacuum is occurring roughly 
instantaneously across space, and in which bubbles of this vacuum collide with roughly equal 
sizes on average. In any different slicing, these vacuum transitions occur at different times across 
space, and bubbles of a given vacuum will collide with systematically different sizes, imparting 
net average momentum to the daughter radiation. In order that the successive generations of 



* Recall that the power spectrum for a field g is defined by (gicffk'} — {2Tr)^6^{k — k')^ Pg 



8 



radiation not encode this momentum, it is then clear that comoving hypersurfaces, on large 
scales, will be slices of constant 0. The coordinate transformation to comoving coordinates, 
with a time shift of St{x.), then induces an additional contribution to the curvature perturbation 
of size A7?.k = H (5tk = Since the fractional change in V during one Hubble time is 

small, this additional term added to TZ in making the transformation is actually much larger 
than the original contribution TZq that we started with. Thus when a mode of the curvature 
perturbation exits the horizon during inflation, it has a size in comoving coordinates given by 

m 



7^k = -ff^. (16) 



Furthermore, it is known that the size of TZ^ in comoving gauge will be time independent 
outside of the horizon [12] . We may thus use expression (TB] to compare our result [12] for the 
power spectrum with the COBE normalization [T5] We obtain 

^'"^ > ~ a{HT) ' 

where H and r are to be evaluated at horizon exit. The following section is thus devoted to 
determining the functions a{HT) and (3{Ht) characterizing the statistical evolution of the field 
at a fixed point in space. 



3 Fluctuations along one Worldline 

Consider a point (x, t) in de Sitter space. We will be interested in this point's entire past 
light cone back to an initial time "t = 0" at which we will say the field value was at step "0" . 
We would like to calculate the statistical properties of the vacuum number occupied by (x, t) . 
Note that our assumption of a uniform initial condition was an arbitrary choice: The random 
fluctuations in the tunneling process will wash out the effects of the initial condition in around 
a Hubble time. What we will be interested in is the asymptotic behavior of (0(x, t)) = {4>{t)) 
and (pRMsi'^^'ty = (pRMsit)'^ for times much larger than 1/H. 

Since we have assumed that the tunneling rate is independent of the step number in the 
chain, the probability per unit volume per unit time for a bubble nucleation event is a fixed 
number, F. The probability of there being N bubble nucleations in the past lightcone of (x, t) 
is then given by the Poisson distribution 

Pr{N,V) = ^{TVre-^'', (18) 
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Figure 2: The past light cone of the point (x, t) in de Sitter space, containing a randomly 
generated set of 7 bubble nucleations. The horizontal axis represents physical distances. 

where V is the 4-volume in the past lightcone back to the initial surface. The mean number of 
bubble nucleations is TV. 

A sample history for (x, t) with bubbles may then be represented by distributing the 
bubbles at random throughout the past lightcone, as shown in figure [2l The question then is, 
given a configuration of nucleation sites in this way, what is the resulting value of the field at 
(x,t)? 

Each nucleation site will have the boundary of its future lightcone traced out by its own 
expanding bubble wall. As shown in figure [3l we could in principle determine the vacuum 
number at (x, t) by following along the sequence of bubble nucleations and wall collisions. This 
picture is rather complicated, and inefficient for the purpose of computer simulations. This 
is essentially because calculating the final field value in this way requires keeping track of the 
vacuum numbers for a continuum of points, and then evolving them with time. In order to 
avoid doing this, we will instead use the following simpler picture: 
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Figure 3: The bubble nucleations from figure [2], shown again with their associated bubble walls 
and vacuum numbers. 



Label each nucleation site "i" with an integer (f){i), corresponding to the vacuum number in 
the interior of the newly formed bubble. Also, let the past lightcone of a spacetime point or a 
nucleation site be denoted by L(x, t) or L{i), respectively. 

Then the vacuum number at (x, t) is given by 



(x,t) = max 0(j), 
ieL(x,i) 



(19) 



where we take the maximum of vacuum numbers over all nucleation sites in the past lightcone 
of (x, t) . Here we have ignored the vanishing probability that there is a bubble nucleation at 
the point (x, t) itself. The vacuum number at an actual nucleation site is given by 



= max (f){j) + 1. 



(20) 



Expression [19] gives us the vacuum number at the chosen point for a particular configuration 
of bubble nucleations. The statistical properties of 0(x, t) will then be determined by summing 
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up over contributions from all such configurations. Obtaining analytic expressions for the 
resulting mean and variance for 0(x, t) is a very hard problem, and for this reason we will 
resort to computer simulations rather than analytic analysis from this point forward. In fact, 
the simple characterization of vacuum numbers given by [19] and [20] immediately suggests a 
useful algorithm for this purpose, which we now describe. In what follows, we will ignore 
certain approximations that were adopted in the code in order to speed up running time. A 
complete discussion of these approximations will be presented in appendix [X] 

To find the vacuum number at (x, t) for a specific history of nucleations, we do the following: 

1. Time order the bubble nucleation sites, and put them in a list. 

2. Associate a vacuum number with each site in the list, initialized to zero. This number will 
be replaced with its correct value as the simulation progresses. 

3. Go through the nucleation sites in order from earliest to latest. At each site, set the 
vacuum number according to equation [201 by searching over all prior sites. 

4. After the vacuum numbers for all sites have been correctly set, use equation [T^] to find the 
vacuum number at (x, t). 

Now, the goal of the simulations is ultimately to calculate (</)(x, t)) and {/)j?a./s(x, t)^. To find 
these expectation values, we have employed a monte carlo method: Dots are distributed at 
random throughout the past light cone of (x, t), in numbers given by the Poisson distribution 
[T51 We then find the resulting </)(x, t) and repeat the process many times, determining (</'(x, t)) 
and (pRMsi'^^ty statistically. 

The results of our simulations are presented in the next section. 

4 Results 

We ran the simulations in 1 + 1, 2+1 and 3+1 dimensions (analytic expressions are trivially 
obtained in 0+1 dimensions), typically for 25k events for each configurationjfl We universally 
find the expected result that both {4>{t)) and (pRusi^Y increase linearly with time. Figures [1] 
through [H] show plots of our data for each dimensionality, in log-log coordinates. Each set of 

^Fractional statistical errors for (0) are much smaller than those for, ^^^ms- therefore also ran smaller 
numbers of events with more calculationally-intensive parameters, which gave good resolution on the former 
but not the latter. For those points, the results for 4>\f^ig are omitted from the analysis. 
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A 


m 


B 


n 


0+1 


1 


-1 


1 


-1 


1+1 


1.45 + 0.06 


-1.001 + 0.002 


0.34 + 0.03 


-0.69 + 0.02 


2+1 


1.55 + 0.35 


-1.00 + 0.01 


0.17 + 0.01 


-0.48 + 0.02 


3+1 


1.58 + 0.30 


-1.00 + 0.02 


0.093 + 0.006 


-0.34 + 0.02 



Table 1: Asymptotic power-law parameters for a and /5, defined in equations [2T] and [221 



measurements of {4>{t)) and (pRMsit^ for the individual r are fit with straight lines (which can 
appear curved in the log- log coordinates of the plots due to nonzero intercept). The error bars 
are typically (6.3 x 10~^)(j)jiMS for {(/)), and (8.9 x 10~^)</)|.^.j-5. for <t>]iMS^ which are both much 
smaller than the data points on the plots. Figures [10] and [11] show plots of the slopes a and (5 
derived from these fits. 

Using equation [TT] we immediately find that nowhere in the range of r's explored by our 
simulation does chain inflation produce an acceptable power of density perturbations, placing 
a modest lower-limit (Hr)^^ ^ 24. In order to estimate the value which reproduces the COBE 
normalization, we perform an extrapolation of our data based on the assumption that, for 
asymptotically large values of {Ht)~^, 

a{HT) ~ AiHr)"" (21) 

(3{Ht) ~ 5(iJr)". (22) 

We find that this ansatz is indeed consistent with the data in all dimensions, and in particular 
in lower dimensions where we ran to larger values of {Ht)~^. Details regarding the extrapo- 
lation method can be found in appendix [A] The results for the parameters characterizing the 
asymptotic behavior are summarized in table [1] Using the above expressions for a and /5 in 
combination with equation [TT] we then obtain the estimate 

(iJr)"^ = (2.2 + 0.8) X 10^ (23) 

This number, when multiplied by the factor A ~ 1.5 gives the required number of vacuum 
transitions per e-folding of inflation. 

Extrapolating beyond the range of parameters directly probed in our simulations was un- 
avoidable, since running the simulations with large numbers of bubble nucleations is extremely 
time and memory consuming. We ran the simulations for several months on a cluster of 32 
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Figure 4: (0) graphed as a function of t in 1+1 dimensions with different r's. Points of each 
r are connected by best-fit straight hues in untransformed coordinates. From bottom to top, 
(Ht)-^ = 3, 6, 12,..., 1536. 




Figure 5: (j)%MS gi'aphed as a function of t in 1+1 dimensions with different r's. Points of each 
r are connected by best-fit straight hues in untransformed coordinates. From bottom to top, 
(Hr)-^ = 3, 6, 12,..., 384. 
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Figure 6: (0) graphed as a function of t in 2+1 dimensions with different r's. Points of each 
r are connected by best-fit straight hues in untransformed coordinates. From bottom to top, 
(Ht)-^ = 3, 6, 12,..., 96. 
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Figure 7: (j)%MS gi'aphed as a function of t in 2+1 dimensions with different r's. Points of each 
r are connected by best-fit straight hues in untransformed coordinates. From bottom to top, 
(Hr)-^ = 3, 6, 12,..., 48. 
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Figure 8: (0) graphed as a function of t in 3+1 dimensions with different r's. Points of each 
r are connected by best-fit straight hues in untransformed coordinates. From bottom to top, 
(Ht)-^ = 3, 6, 12,..., 48. 




Figure 9: (j)%MS gi'aphed as a function of t in 3+1 dimensions with different r's. Points of each 
r are connected by best-fit straight hues in untransformed coordinates. From bottom to top, 
(Hr)-^ = 3,6,12,24. 
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Figure 10: The function a, as described in the text, plotted versus {Ht)~^ for all 3 dimensions. 
For clarity, the points for 1+1 dimensions and 2+1 dimensions have been offset vertically by 
factors of 25 and 5, respectively. 
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Figure 11: The function /3, as described in the text, plotted versus (Ht) ^ for all 3 dimensions. 
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Figure 12: Spatial 2-point correlator of as a function of separation. This was generated by 
the simulation with 1+1 dimensions, (ifr)~^ = 6, Ht = 16. The slope of the best-fit line is 
-1.259 ± 0.022. The /3 value for this set of parameters is 1.262 ± 0.006. 



IGHz processors. However, it is clear from the analysis of appendix lAl t hat for any given dimen- 
sionality, a asymptotes to ~ {Ht)~^, and /3 is consistent with power-law behavior. Of course 
it is possible that a drastic change in the behavior of the system could take place beyond the 
range of our simulations, but given the simple scaling evident in figures [TU] and [TTl this seems 
rather unlikely. 

Finally, for a few isolated points in parameter spacej^ we ran full simulations for the 2-point 
correlator (0(x, t) 0(y,t)) — in order to show agreement with our expression [TO] for how 

it should be determined in terms of (pmisi^Y ■ A typical result is shown in figure [T2j Statistical 
agreement with the expected slope was found for all parameter points that we checked. 

^Specifically, in 1+1 we checked (Ht)^^ = 6 and 8, Ht = 8 and 16, all four combinations. In 2+1 we checked 
(Ht)-^ =3, Ht== 16. 
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5 Conclusions and Discussion 



We have demonstrated that Chain Inflation, with the assumption of uniformity of tunnehng 
events at each step of the chain, requires about 10^ vacua per e-folding of inflation. This is 
necessary in order to produce density perturbations of the size measured by COBE. 

Perhaps the primary model in which one might imagine having a uniform sequence of 
tunneling events is the axion model of [7]. The axion potential takes the form of a tilted 
cosine, with the number of local minima scaling like the number of quark flavors charged under 
the Peccei Quinn symmetry [13]. We have thus shown that this model produces unacceptable 
density perturbations, unless the number of fermions in the theory is very large. 

While our result for the density perturbations does depend on the assumption of uniformity 
of tunneling events, we may still put a constraint on more general setups for chain inflation. It 
is clear that for a flxed total number of transitions, the effect of allowing the tunneling to be 
non-uniform will be an increase in the size of the perturbations, not a decrease. Thus in the 
non-uniform case, 10^ vacua per e-folding of inflation becomes a lower bound. 

The large number of vacua we have shown to be required by the model is worrisome, as 
it seems hard to imagine such a long sequence of vacua being able to exhibit the properties 
needed for successful chain inflation. For example, it is necessary that the system never enter 
a stage of eternal inflation. Moreover, since about one e-folding worth of radiation is left over 
from previous percolations at any given time, the energy in the radiation bath will be about 
10^ times the vacuum energy difference between successive steps of the chain. The effects of 
such a large quantity of radiation on the evolution of the system are thus very important to 
explore. 

There are clearly many open problems in chain inflation, and much room for future work. 
Our results were obtained through extrapolations of computer simulations, and trying to flnd 
an analytic understanding of the power law behaviors [2T] and [22] seems important, although 
the problem has proven rather difficult. It is also an open question as to whether or not 
chain inflation could yield observable non-gaussianities in the CMB. Finally, there is another 
possibility worth considering: It might seem that the assumption of uniformity of tunneling 
events was necessary in order to obtain a scale invariant spectrum of fluctuations. On the other 
hand, since we know that many vacuum transitions are required per e-folding of inflation, a 
given measurement of a CMB multipole moment actually averages over the results of a large 
number of tunneling events. It is conceivable that even if the tunneling rates were varying, they 
might be sufficiently uniform "on average" to produce a scale invariant spectrum. It would be 
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interesting to explore this possibility, perhaps with a string-inspired chain of much more than 
10"^ vacua. 
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A Computer Simulation 



Here we shall discuss in more detail the algorithm used in our monte carlo simulations. 

The volume of a simulation is a l)-dimensional cube of physical side length 2H^^ crossed 
with an interval of length t, the time extent of the simulation. This is meant to represent a 
D-volume in unperturbed de Sitter space in physical coordinates, essentially a region of fixed 
spatial extent about a particular free-falling worldline centered at the origin. The lightcone of a 
point in this coordinate system can be time-sliced into (D— l)-spheres of radius \e^^^ — l] , 
centered about the spatial position of a co moving worldline passing through that point. The 
set of points with |x| = specifies the causal horizon. For any point within this spacetime 
"cylinder," the past lightcone will be completely contained as well. 

We construct this volume to investigate the field fluctuations at the spacetime point (0,t). 
The D-volume described above contains the entirety of this point's past lightcone in that time 
interval, which is restricted to lie within the horizon. Assuming that inflation starts at t = 0, 
this D-volume then contains all bubble nucleations that will determine the field value, as well 
as their past lightcones. An illustration of the shape of (0,t)'s past lightcone can be seen in 
figures [2] and [31 

Since the bubble nucleation density per unit /^-volume is constant for all tunnelings, we 
distribute nucleation sites in a statistically uniform manner across the volume with no prior 
knowledge of which vacuum is nucleating at each site. First, we throw the number of sites 
from a Poisson distribution with mean Vt {2H~^)^ ^. We then place these sites at uniformly 
random, uncorrelated positions in the D- volume, disregarding those nucleations that fall outside 
of (0,t)'s past lightcone. 

To calculate the statistics of </)(0,t), we implement the procedure described in section [3] 
for a large number (typically tens of thousands) of random nucleation configurations. Though 
following this basic algorithm, the simulation also incorporates a few time-saving tricks. 

First, the past lightcone search for each nucleation site, necessary to implement equation [2Ul 
is cut off at a time tcut earlier. This is motivated by the fact that successive percolations "screen" 
each other. The parent bubble in which a given nucleation will take place will typically have 
been birthed no more than a few r in the past. It is obviously counterproductive to consider 
nucleations that occurred many r earlier, as the bubbles generated during those early times 
will have been crowded out by more advanced bubbles. Consequently, we set tcut to be a few 
r, typically 3r for D = A, and 5r for = 3 and D = 2j^ We find that these values yield 

^In lower dimensions, a point in spacetime sees a smaller number of nucleations in its recent past for fixed 
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statistical results that are stable under tcut increases. 

The second trick is to divide (D— l)-space into numbered cubical "boxes" in order to organize 
the nucleation sites. The boxes are O(tcut) along each side, so that we can find all sites in a 
past lightcone out to tcut by searching just within boxes in the immediate vicinity. Distant 
sites are thus automatically ignored. More rigorously, we calculate the spatial extent of the 
past lightcone tcut back, and search over all boxes that might contain it. Since this procedure 
captures the entire past lightcone out to tcut^ it does not constitute an additional approximation. 

We also perform a handful of modified simulations in order to calculate the 2-point correlator 
{(f){A) (j){B)) for spacetime points A and B outside of each other's respective horizons. These 
follow largely the same procedure, with simulations for A and B run in parallel. Points that 
actually share history will have common nucleations in their past lightcones. Consequently, 
some of the nucleation sites in A's and S's simulation volumes are actually "the same," and 
we correlate their placement accordingly. 

B Extrapolation Method 

Here we present the analysis techniques used to extract the results for table [T] in section |H 

We ran the monte carlo simulation in D = 2, 3, 4, with nucleation rates {Ht)~^ = 3, 6, 12, . . ., 
and for time durations Ht = 1, 2, 4, 8, . . ., as far as our processor and memory resources would 
allow within a reasonable time. The set of measurements for {4>{t)) and (pRMsitY appear in 
figures m through [91 We fit these to obtain the slopes a^Hr) and i3{Ht), presented in figures 

In D = A, the largest {Ht)~^ for which our statistics allowed both a and (3 to be precisely 
calculated was 24. Equation [T7] indicates that the resulting values are far from reproducing the 
correct scale of density perturbations. Consequently, we are forced to perform an extrapolation 
to much larger values of (Hr)^^. We assume that a and /5 asymptote to the power law forms 
of equations [21] and [22l In order to estimate the asymptotic power laws, we first fit ascending 
pairs of (Hr)^^ points to power laws, and then investigate the convergence of the coefficients 
and powers returned by the fits. The results of these pairwise fits are illustrated in Figures 
[T^ through [ini Note that for a, points with different D's are vertically offset for clarity of 
presentation. In reality, they essentially lie on top of each other. 

For a, the statistical errors are small and the convergence behavior is clear. However, for 

tcut than a point in higher dimensions would, and thus has a greater chance of missing the nucleation(s) with 
maximum (f>. 
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Figure 13: Best-fit power-law coefficients for a using pairs of points {Ht)~^ and 2{Ht)~^. For 
clarity, the points for 1-1-1 dimensions and 2-1-1 dimensions have been offset vertically by 0.2 
and 0.1, respectively. The extrapolations (dashed lines) are described in the text. 
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Figure 14: Best-fit power-law powers for a using pairs of points {Ht)~^ and 2{Ht)~^. For 
clarity, the points for 1-1-1 dimensions and 2-1-1 dimensions have been offset vertically by 0.04 
and 0.02, respectively. The extrapolations (dashed lines) are described in the text. 
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Figure 15: Best-fit power-law coefficients for (3 using pairs of points {Ht) ^ and 2{Ht) ^. 
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the physically relevant D = 4 case, we were not able to go far enough in {Ht)~^ to explicitly see 
the convergence as well as we can in lower dimensions. In order to perform the extrapolation, 
we fit the power-law parameters in all 3 D's with functions of the form pi + p2{HTy^ , where 
the Pi are the fit parameters. These are seen to be visually good matches in the figures, though 
note that we do not claim that they furnish accurate models of the convergence behavior. 
(The /NDF values for these fits range from 1 to 100.) The best-fit asymptotic values pi 
constitute our central estimates. We then define the systematic error on each of these as the 
difference between the extrapolation estimate and the power-law parameter measured for the 
last pair. Our estimate for the D = A asymptotic coefficient and power are thus 1.58 ± 0.30 
and —1.00 ± 0.02, respectively. Interestingly, the results from all dimensions taken together 
suggest a universal asymptotic form a(iJr) ~ (1.5) (i/r)^^, though we do not assume this in 
our analysis. In this limit, {4>(t)) is independent of H. In any event, the actual dependence on 
H is seen to be very weak. 

For (3, the statistical errors are too large to resolve any details of the hypothesized conver- 
gence behavior. All pairs of points yield consistent coefficients and powers for the power-law 
fits§ We therefore conservatively define the central estimates and their errors as the parameter 
values and errors obtained from the power-law fit of the last pair. For D = 4, the estimated 
asymptotic coefficient and power are 0.093 ± 0.006 and —0.34 ± 0.02, respectively. 

Another source of potential systematic error comes from the tcut cutoff on past lightcone 
searches, described in appendix |X1 We examined the numerical stability of our a and (3 values 
with respect to variations in tcut in all D for several values of r, and found convergence to 
within the statistical errors. Thus, we do not associate an error with this variable. If it is there, 
it is sub dominant. 

The remainder of our estimated asymptotic power-law forms for a and /3 are summarized 
in table [H 



*The D = 3 points do display some slope by eye, and a fit to a flat line indicates that convergence across the 
entire sampled range of {Ht)~^ is only about 1% probable. However, consistency with a flat line does emerge 
past the first point. 
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